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The effective interaction between a sphere with an open cavity (lock) and a spherical macroparticle 
(key), both immersed in a hard sphere fluid, is studied by means of Monte Carlo simulations. As 
a result, a 2d map of the key- lock effective interaction potential is constructed, which leads to the 
proposal of a self- assembling mechanism: there exists trajectories through which the key-lock pair 
could assemble avoiding trespassing potential barriers. Hence, solely the entropic contribution can 
induce their self-assembling even in the absence of attractive forces. This study points out the 
solvent contribution within the underlying mechanisms of substrate-protein assembly /disassembly 
processes, which are important steps of the enzyme catalysis and protein mediated transport. 



The key-lock (KL) self- assembling mechanism [l| is 
found in several vital biological processes such as ac- 
tive protein mediated transport Q , enzyme catalysis Q , 
DNA/RNA transduction and replication among oth- 
ers. It consists on the perfect match of a macromolecule 
referred to as the key, into an irregular open cavity of 
another generally larger macromolecule, i.e., the lock. In 
case of catalysis, the lock particle is an enzyme, which 
may have more than one cavity to capture different reac- 
tants (substrates). The perfect match between the sub- 
strates and the enzyme guaranties, at least partially, the 
specificity of the desired reaction, since the active site 
(the catalyst active part of the enzyme) is generally sit- 
uated inside the lock. The specificity is so high that 
scientists have tried to emulate it for designing their own 
catalysts [5KZj. 

In general, the catalatic and protein mediated trans- 
port kinetics involve the following steps: assembling, re- 
action or transport, and disassembling [8]. Since these 
steps are sequential, a slow step would strongly affect 
the overall rate of the whole process (bottleneck). Most 
biological catalytic and protein mediated transport pro- 
cesses are not only specific but also show huge kinetic 
rates 0, strongly suggesting that the complicated 
match between the key and the lock and its dissembling 
do not retard the kinetics of the whole process. 

This communication focuses on the collective contribu- 
tion of the solvent to the effective interaction potential 
between the key and lock macromolecules, lea ving aside 
the evidently important role of the electrostatic [lOj , Lon- 
don - van der Waals, hydrogen bonding, end internal 
KL entropy contributions [11[. Additionally, the consid- 
ered KL / solvent size ratio and solvent density are still far 
from a realistic biological system, and so, in this context 
the results reported here must be considered as a guide 
only. Since just the exclusion potential energies of the 
solvent particles and macroparticles are accounted for, 
the resulting macroparticles interaction force is referred 



to as entropic [Hj, depletion [l3|, or contact 15, 
Although solely this interaction is capable of producing 
phase transitions, clusters growing, and self-assembling 




FIG. 1. a) Schematic representation of the studied system 
(y =0 plane cut), b) Snapshot of an equilibrium configura- 
tion; look particle at the left (blue), key particle at the right 
(red), and solvent particles (semitransparent grey spheres). 



[l3L flil . [l7L flil , it is generally ignored by studies of cat- 
alytic and protein mediated transport processes. 

Previously, the pair potential between hard macropar- 
ticles immersed in a hard sphere fluid has been studied 
[T9M23II . In most cases the interaction between convex 
or planar surfaces has been addressed, and only a few 
works deal with concave surfaces [14, 21, \2m. In this 



work we gain insight into the KL self- assembling mech- 
anism by studying the pair potential between a spheri- 
cal macroparticle with an open spherical cavity (concave 
surface), the lock, and a hard sphere that fits the cavity, 
the key, both immersed in a bath of smaller hard sphere 
particles. To achieve this goal it becomes necessary to 
analyze the whole 2d energy map from where low energy 
trajectories can be deduced. These trajectories, which 
may avoid trespassing potential barriers, are those which 
would make possible a fast KL self- assembling kinetics. 
As mentioned, this is a characteristic of the enzymatic 
catalysis and the protein mediated transport. 

Three species are considered: the lock, consisting on 
a hard sphere of radius n=3cr having a spherical cavity 



of radius rh located at a distance d = \f~r 



from 



2 




; 



0.5 
0(TT) 



1.0 



FIG. 2. Solvent density profile at the key particle surface, p sc , 
as a function of the angles and <f> for (x,z) = (0.2<j, 1.6a). 
The coordinate system {x',y',z') is located at the center of 
the key particle. Darker colors mean a higher value of p sc ] 
white means zero. 



its center; the key, a hard sphere of radius rk<rh] and 
the solvent, a hard sphere fluid made of particles of di- 
ameter cr<r/c. Thus, the KL closest approach distance is 
dki = d — rh + rk- The center of the lock particle is fixed 
at (x,y,z) = (0,0, — dfcz)> being the origin of coordinates 
placed at the center of the simulation box. The position 
of the key particle is also at the y = plane, at a given 
(x, z) point which varies from run to run. Particles and 
locations are schematized in Fig. [TJi. Unless otherwise 
indicated, computer experiments were done by setting 
^=^fc=l-7cr, so that dki=d. The lengths of the simu- 
lation box sides are L z =18a and L x =L y =12a. Solvent 
particles are initially randomly placed, and then moved 
according to the Monte Carlo scheme. The bulk solvent 
volume fraction is set to p so =0.2. A snapshot of the sys- 
tem is shown in Fig. 

The contact force acting on a macromolecule is given 
by F=—kBT f A p sc iids, where p sc is the density profile 
of the solvent at the macromolecule surface, n is a unit 
vector pointing out the surface, ks is the Boltzmann con- 
stant, T is the temperature, and the integral subindex A 
refers to the macroparticle surface area. p sc for the key 
particle located at (x, z) = (0.2<j, 1.6<r) is given in Fig. [2] 
as a function of the angles and <j). These angles are 
defined as shown in the figure. For the surface region far 
from the lock, i.e., > 0.57T, p sc is practically constant, 
p sc « 3.4p so (lighter yellow). Only the region close to 
the lock particle shows clear variations of p sc . There is a 
depletion region (white band at « 0.257r) close to the 
cavity border of the lock, where solvent particles cannot 
access. As can be seen, this region (and the whole sur- 
face) is <j) dependent since the axial symmetry around 
z was broken by setting x^O. For < 0.27T, p sc pro- 
duces larger values than those for > 0.57T, since the 
solvent particles are strongly adsorbed into the lock cav- 
ity. Moreover, the largest values of p sc ~ 11.2p so (darker 
red) correspond to particles inside the cavity and close 
to its border. Hence, the region where the cavity and 
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FIG. 3. a) F z as a function of z for x=0. The dotted line 
represents the UHS contact. Solvent densities at the key par- 
ticle surface (as given in Fig. [2]) are also shown for the cases 
indicated by arrows, b) Effective pair potential interaction 
energy, E. c) Closeup of the same energy data. In all plots, 
square symbols are for the KL pair while circles are for the 
UHS pair. 



key surfaces are close is preferred by the solvent parti- 
cles. This suggests that the system is optimizing space 
by placing solvent particles where a large surface / volume 
ratio is found. For large z, p sc tends to be independent 
of and 0, implying that F asymptotically decays to 
zero [25]. This long range behavior of the contact force 
may explain how certain enzymes act on substrates at 
kinetic rates that approach the encounter rate of the KL 
in solution [9j. 

For the case of Fig. [2j the depletion region does not 
contribute to counterbalance the force exerted by the 
solvent on the outer surface of the key, which produces 
a negative (leftwards/attractive) contribution to the z 
component of the force, F z . On the other hand, the fringe 
at < 0.27T tends to counterbalance (and may even over- 
balance since p sc is large in this region) the leftwards F z 
contribution. Thus, in this case the sign of F z is not evi- 
dent, whereas F y =0 by symmetry, and F x ^ 0. Note that 
the average force on the key particle must be equal (with 
opposite sign) to that acting on the lock particle. This 
was verified to guaranty the correctness of the algorithm. 

F z calculated for x=0 is shown in Fig. [3^ as a func- 
tion z (errors are always smaller than symbols). The 
F z which results from substituting the lock particle by 
a hard sphere of radius r\ is also included, i.e., for the 
lock particle without the cavity (unsymmetrical hard 
spheres, UHS). In Figs. [3b and c it is shown Eq = 
(— F z dz)\ x= o, i.e., the corresponding effective pair po- 
tential energy. For both cases F z is attractive at contact. 
However, the cavity induces a KL attractive force more 
than twice larger than that for the UHS. This result is in 
agreement with the theoretical prediction of Kinoshita et 
al. [14|, |2l|. In terms of energy, this turns into a decrease 
of ~14/c#T, i.e., a potential well ~7 times deeper. Hence, 



if the UHS potential well is capable of segregating large 
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particles from the solvent, as found in many colloidal 
systems [26j, then the KL potential should produce its 
self-assembling. 

Forces can be understood by analyzing the solvent den- 
sity profiles at the key surface, p sc . The insets of Fig. [3^ 
compare the p sc obtained for a UHS pair with those ob- 
tained for the KL pair. By symmetry, these panels are 
independent of (p. For the panels corresponding to the 
macroparticles closest approach distances (1 and 2), it 
is clearly seen that the depletion region becomes much 
larger for the KL pair (white fringe), explaining the much 
larger attraction found for it. On the other hand, the 
other two panels are obtained for z=2.6a (3), KL pair 
case, and a surface-surface separation of 1.4a, UHS pair 
(4). For these separation distances the net forces are 
positive (repulsive), producing potential barriers. It is 
observed a single fringe of large p sc values for the UHS 
case, while a double fringe structure appears for the KL 
pair. This last structure is produced by the solvent par- 
ticles adsorbed near the cavity border, inside and outside 
it. In both cases (3 and 4), the fringes overcompensate 
the outer solvent pressure yielding a net repulsive force. 
Thus, the double fringe structure leads to a larger repul- 
sive force than the single one (see the higher peak force 
for the KL pair than for the UHS one in Fig. [3^). 

Both pair potentials show energetic barriers (Figs. [3b 
and c). The KL pair potential peak is at z~2.6a, while 
for the UHS is at 2^3. 6a, i.e., for a surface-surface sep- 
aration distance of «1.2cr (see Fig. [3b). Although the 
peaks are not very high, less than ksT, they may hin- 
der the self-assembling kinetics. It should be pointed out 
that the KL potential barrier is almost twice higher than 
the UHS one. In addition, these barriers enlarge by in- 
creasing the solvent density, so that, they may be quite 
large for denser solvents such as water. Based solely in 
these results, one could conclude that the self-assembling 
kinetics for the KL pair is relatively slow (at least slower 
than the UHS assembly). As we will show, this is not 
true. 

By computing E(x,z) — — f^F z dz'\ x along different 
values of x = const trajectories, the pair potential energy 
for the KL pair as a function of x and z can be obtained. 
The results were fitted by a 2d-function by means of the 
Levenberg-Marquardt method. This function is plotted 
in Fig. HI As expected E(x, z) is not radially symmetric. 
This implies the presence of KL noncentral forces, which 
translates into opposite torques acting on the key and the 
lock macroparticles [12| . For x = torques disappear and 
the trajectory along the map corresponds to the squares 
curve of Fig. [3b. The most important finding is that the 
potential barrier disappears for z trajectories passing at 
\x\ ~ 1.3a. For ip > 7r/10 (see Fig. QJi for the definition 
of ip) the barrier reapers and the energy map converges 
to a radial dependent map as that corresponding to the 
UHS case for ip > 7r/4 (circle symbols of Figs. [3b and c). 
This means that the UHS map is equal to that for the 




FIG. 4. Pair potential energy, E(x, z), for the KL pair. Bluish 
colors mean attraction and reddish colors mean repulsion. 
Darker tones denote larger energy values. Dotted and dashed- 
dotted lines represent low energy trajectories. 
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FIG. 5. a) F z as a function of z for x—0. Effective pair 
potential interaction energy, E. c) Closeup of the same energy 
data. Square symbols represent rh=rk = l>7a, circles rh=l>7a 
and rfc=1.0<7, and triangles r^=2.4cr and r^ — \.7G. 



KL pair with the cavity pointing backwards. 

Low energy trajectories, which do not trespass poten- 
tial barriers, are those shown as dotted lines (due to the 
z axis symmetry these trajectories form a surface of rev- 
olution in space). These trajectories avoid trapping sol- 
vent particles inside the cavity, which would exert pos- 
itive pressure on the KL pair, providing them a way to 
escape. Upon the key particle overcomes the central po- 
tential barrier through the sides, it gets in contact with 
the border of the lock cavity. Then the macroparticles 
slides on the cavity surface towards the (x, z) — (0, 0) 
position. Another relatively low energy way towards this 
position would be trespassing the potential barrier far 
from the cavity and follow the dotted-dashed line trajec- 
tories shown in Fig. [4j As shown by Figs. [3b and c, the 
potential barrier for large tp is smaller than that for small 
ip. Although this way the KL pair still needs to trespass 
a potential barrier, it is likely to occur since the external 
surface of the lock is larger than the cavity area. 
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Up to here the results show that a KL perfect match 
(Th=rk><r) not only produces a very low pair potential 
energy at the closest approach distance but also yields 
low energetic trajectories. These two facts guaranties a 
fast KL self assembling. Hereinafter, we focus on how an 
imperfect match changes the KL pair potential. Two 
other cases are studied; these are: (a) r^=1.7cr and 
rk =1.0a, i.e., the key size was decreased and the lock 
particle is kept as it was, and (b) r^=2.4cr and 7^=1. 7a, 
i.e., keeping the same key the lock cavity is enlarged (see 
the drawings of Fig. [5j). Case (a) may correspond to the 
decrease of the key after a breakup reaction inside the 
lock (an enzyme catalyzed reaction), whereas case (b) 
may model a protein conformation change (for instance, 
after the transportation of the key). As mentioned above, 
the sequential mechanism ends up with the key release, 
so that the catalyst or the transporter can be recovered. 
A priory this seems unlikely due to the large energy well 
found for the perfect match case. 

Results for cases (a) and (b) are shown as circles and 
triangles in Fig. [5j respectively, as well as the KL perfect 
match results (squares). As can be seen, cases (a) and 
(b) reach pair potential wells which are approximately 
half deep than the perfect match case. Hence, both situ- 
ations, a size decrease of the key by following a breakup 
reaction or the enlargement of the cavity produced by a 
conformation change of the lock macroparticle favor the 
key release. 

In summary, our results show that the KL self- 
assembling is favored by the entropic solvent contribu- 
tion. In addition, low energy KL self assembling trajecto- 
ries which avoid trespassing energy barriers are obtained, 
guarantying a fast KL self-assembling kinetics. During 
the process, a net torque appears acting on both parti- 
cles and guiding them to the match position. Finally, a 
lock conformation change (frequently found during catal- 
ysis and active transport processes [U() produces solvent 
entropic contributions favoring the key release. Thus, the 
results agree with the fast asembling/dissasembling steps 
which occur during the protein mediated transport and 
the enzyme catalysis processes. Finally, we emphasize 
that works on enzyme and transport kinetics should ac- 
count for the solvent entropy changes. 
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